Phase transitions in nanosystems caused by interface motion: The Ising bi-pyramid 

with competing surface fields 

A. Milchev,i'2 M. Muller,!^^ ^j-^j Binder^ 

^ Institut fiir Physik, WA 331, Johannes Gutenberg Universitdt, D 55099 Mainz, Germany 
^Institute for Physical Ghemistry, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria 
^Department of Physics, University of Wisconsin-Madison, 
1150 University Avenue, Madison, WI 53706-1390 

The phase behavior of a large but finite Ising ferromagnet in the presence of competing surface 
magnetic fields zb/Zj, is studied by Monte Carlo simulations and by phenomenological theory. Specif- 
ically, the geometry of a double pyramid of height 2L is considered, such that the surface field is 
positive on the four upper triangular surfaces of the bi-pyramid and negative on the lower ones. It is 
shown that the total spontaneous magnetization vanishes (for L —> oo) at the temperature Tf{H), 
related to the "filling transition" of a semi-infinite pyramid, which can be well below the critical 
temperature of the bulk. The discontinuous vanishing of the magnetization is accompanied by a 
susceptibility that diverges with a Curie- Weiss power law, when the transition is approached from 
either side. A Landau theory with size-dependent critical amplitudes is proposed to explain these 
observations, and confirmed by finite size scaling analysis of the simulation results. The extension 
of these results to other nanosystems (gas-liquid systems, binary mixtures, etc.) is briefly discussed. 

PACS numbers; 68.08.Bc, 05.70.Fh,68.35.Rh, 64.60.Pr 

I. INTRODUCTION 

The current paradigm of attempting to develop various kinds of nanoscopic devices requires careful consideration 
of the phase behavior of nanosystems, since in nanoscopic geometries effects due to external walls or other boundaries 
of the system can modify its "bulk" behavior substantially. Qualitatively new kinds of phenomena may occur, that 
have not yet been studied for macroscopic bulk systems. 

We demonstrate a new kind of phase transition in the present paper, which belongs to the class of interface 
localization-delocalization phenomena, using the simple Ising ferromagnet with nearest-neighbor exchange on a cubic 
lattice as a generic example. Choosing a compact octahedral shape of the system in the form of a bi-pyramid of height 
2L, we assume that on the upper surfaces of the pyramid {0 < z < L) a positive surface magnetic field -\-Hs acts, while 
on the lower surfaces (with —L < z < 0) the field is negative but of the same absolute strength, so that no sign of the 
magnetization is overall preferred. More generally, one might consider the case with positive and negative fields of 
different strength; their difference, however, could be effectively compensated by a suitably chosen bulk field such that 
at low temperatures again a degeneracy with respect to the sign of the spontaneous magnetization is possible, similar 
to the case of "capillary condensation"- type phenomena in semi- infinite thin films 1, 2, 3, 4, 5, 6,, 0, |M S Eo, H]- In 
this case one can also expect an interesting interplay between the wetting behavior of the semi-infinite system and 
the phase behavior in confinement, a complication that is not considered in the present manuscript. 

Such a system is then described (for L oo) by an order parameter (the spontaneous magnetization of the 
Ising ferromagnet), which does not remain non-zero up to the critical temperature Tct of the bulk three-dimensional 
model, but rather only up to a temperature Tf( Hs), identical with the (critical) temperature of the filling transition 
I3IIllIlliilIlllElIl|IllM|2l|2^ ElM IM Hi in a single semi-infinite pyramid. As will be 

discussed in this paper, this new kind of phase transition [3l| in the limit L —> oo can be either of first order or of 
second order, depending on the value of the line tension, which describes |33 . l33.l33 .r35. 36. 37] the free energy excess 
associated with the contact line where the interface separating oppositely oriented domains meets the free surface (or 
inert wall that confines the system, respectively) . Of course, as long as the linear dimension L of the system is large 
but finite, finite-size rounding of this phase transition needs to be considered, and hence we shall present a tentative 
generaHzation of the theory of finite size scafing |33, ISa, EO) EJ) 113) EM IMi 113 to the present situation here. 

A qualitative explanation of this new transition is sketched in Fig. ^ It is assumed that the surface magnetic field 
strength Hg is small enough, so that for zero temperature the ground state of the system has a uniform (positive or 
negative) magnetization, in spite of the unfavorable energy cost due to the surfaces (or walls, respectively) where the 
surface magnetic field is oppositely oriented to the direction of the spontaneous magnetization. As the temperature is 
raised, the interfacial free energy a between oppositely oriented domains decreases faster than the excess free energy 
difference, fs {Hg , T) , of a positively oriented domain between surfaces with ±Hs . As is well known the interfacial 
free energy a vanishes at the bulk critical temperature Tct according to a power law cr (x (1 — T/TcbY'^'' with the 
correlation length critical exponent 0,^3 ^fe ~ 0.63. At the temperature Tf{Hs), these surface free energies become 
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FIG. 1: Ising ferromagnet on a simple cubic lattice whose surfaces form a bi-pyramid (left) and resulting phase transition in the 
limit L oo (right), plotting the spontaneous magnetization ms{T) versus temperature T. Signs (+, — ) along the cross section 
of the bi-pyramid (left upper part) or on the triangular projections of the surfaces in the top view (left lower part) refer to 
the surface magnetic field, ±Hs, that acts on the spins in the surface planes only. Note that the basal plane of the bi-pyramid 
is taken to be the (xy) plane of the simple cubic lattice, and measuring lengths in units of the lattice spacing, each pyramid 
takes L planes (with a single spin in the pyramid top), so the total linear dimension from top to bottom of the bi-pyramid is 
2L (the extra lattice slice accounts for the basal plane common to both pyramids). For T < Tf{Hs), the interface between 
the domains with negative (j.) and positive (t) magnetization is located close to one of the corners (e.g. the bottom corner, as 
assumed in the figure, see the magnified view). The local fiuctuations of the interface extend over a correlation range as 
indicated by the double arrow. As the temperature is increased towards the filling transition temperatures, T — > Tf{Hs), the 
interface detaches from the corner and moves towards the midplane of the bi-pyramid. For T > Tf[Hs) the magnetization rris 
then remains zero. 

equal, and hence for T > Tf{Hs) it is energetically favorable to have a state with two oppositely magnetized domains, 
separated by an interface located in the basal plane of the bi-pyramid (Fig. Of course, in the actual calculation 
of Tf{Hs) not only the surface free energies tr and fs{Hs,T) (per unit area) matter, but one must also consider the 
fact that the four triangular surfaces take a larger area (depending on the opening angle a of the pyramids) than the 
area of the interface. In this work, we consider explicitly only the case a = 45°, so that the surfaces of the bi-pyramid 
meet at the basal plane at an angle of 90° , but it is clear that the general features of the phenomena described here 
do not depend on this particular choice. In fact, we speculate that also the choice of planar surfaces is an irrelevant 
detail, and similar behavior could be observed for other geometries such as double-cones with different surface fields 
on the upper and lower portion. 

The outline of this paper is as follows. In Sec. II, we recall the basic facts about the filling transition in semi-infinite 
cones, and develop a tentative phenomenological theory for describing the transition explained in Fig. ^ Sec. Ill 
describes our Monte Carlo results and interprets them in terms of the phenomenological description of Sec. II. Finally 
Sec. IV contains our conclusions, and discusses briefly the extension to phase transitions of other systems (gas-liquid 
systems, binary mixtures, etc.) in related geometries. 

II. THEORETICAL BACKGROUND 

In this section two complementary phenomenological approaches to the phase transition in a double-pyramid are 
developed. In subsection A we use the description of the filling of a single cone '2l'|. This yields the location of the 
filling transition in the limit L ^ og and we discuss modifications due to the double-pyramid geometry. This approach 
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FIG. 2: Comparison of a bi-pyramid of height 2L and basal plane of linear dimension 2L with a single (infinitely large) 
pyramid, assuming the same positive surface field Hs at the corresponding surfaces, and considering the situation that the 
interface between the domains with positive (}) and negative {[) magnetization is located at the same distance io from the 
bottom corner in both cases. An opening angle a — 45" is assumed for simplicity. 



is expected to yield a good description if the magnetization is close to its saturation value, i.e., for it ^ 1, where 
t denotes the reduced distance from the filling transition. The role of fluctuations within this context is considered 
in subsection B. Then, in subsection C, we develop a phenomenological Landau-type theory for the case that the 
interface fluctuates around the basal plane. This approach is able to describe the behavior in the ultimate vicinity of 
the transition, L'^t <C 1, and the fluctuations above the transition. 



A. Phenomenological considerations in terms of surface thermodynamics 

Our phenomenological description assumes that the theory of cone filling can be directly applied to the filling 
of a semi- infinite pyramid (i.e., we ignore the excess free energy at the edges of the pyramid, where the contact lines 
of the interface with two triangular pyramid surfaces meet). We compare the bi-pyramid geometry to an equivalent 
situation of a semi-infinite single pyramid, and consider the case when the interface is located at a height €o above the 
bottom corner (Fig.|2Jl. We write the free energy of the semi-infinite pyramid, relative to a state with no interface, in 
terms of surface and line free energies 

A^, = Alia + 8£o^Tiinc - A^/2^lf,{H,) (1) 

In Eq. ^ we have used the geometrical factors appropriate for the opening angle a = 45° (a generalization in terms 
of other choices for the angle a is straightforward), and we have also suppressed temperature arguments throughout, 
in order to simplify the notation. Actually, rather than using the temperature, T, as a control parameter as assumed 
in Fig. n] we find it more convenient to use the strength of the surface magnetic field Hs instead. (In the plane of 
variables T, Hg the filling transition line is described by the inverse function HsciT) of the function Tf{Hs). As long 
as one crosses this line under a finite angle, it does not matter whether T or Hs is used as a control variable). 

Since we know that at the filling transition the interface can move infinitely far apart from the lower corner, Iq oo, 
we must have AF = for Hg = Hgc, i.e. 

a = V2f4Hgc). (2) 

This result agrees with the macroscopic filling condition that the cone fills if the contact angle on a planar substrate 
equals the cone angle, a. In the vicinity of Hsc the variation of fs{Hs) with Hs is linear, fs{Hs) — fs{Hsc) + {Hs — 
Hsc)fs- As a result, near Hs = Hsc Eq. can be rewritten as 

AF,(4) = 84aii„c - AV2el{Hs - Hsc)f's (3) 
Minimization of Eq. ^ with respect to £q readily yields 

= '^""'^ (4) 

V2{Hs - Hsc)f's 

Since fs{Hs) is a monotonously increasing function of Hs,fs > and hence the denominator of Eq. Q is negative 
in the considered region Hs < Hsc (for Hs > Hsc the pyramid is "filled", i.e., £o = on this other side of the filling 
transition). Of course, only non- negative solutions for ig are physically meaningful, and hence we require that the 
line tension is negative, crnne < 0. The fact that critical cone filling can only occur for negative values of the line 
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FIG. 3; Schematic plot of the surface free energies of a large but finite bi-pyramid Ising system versus Hs (or t = Hs/ Hsc — 1, 
respectively). For L — > oo the surface free energy is since the contribution due to the interface 

at height £o above the lower corner is negligible. For Hs near Hsc the variation of fs{H, T) with Hs is linear. For Hs > Hsc the 
surface free energy is due to the interface in the basal plane of the bi-pyramid, 4a L^, independent of Hg, in the limit L ^ oo. 
If large but finite linear dimensions L are considered, the surface free energies are reduced because of line tension contributions. 
For Hs > Hsc(L), which is characterized hy Iq — L (dot in the figure), this reduction is 8a{i^^L (in the vicinity of Hsc the 
dependence of cr{i^^ on Hs can be neglected). For Hs < Hsc{L), the depression AFs of the surface free energy relative to its 
asymptotic expression for L ^ oo gradually grows with increasing Hs, reflecting the gradual increase of to. Assuming that the 
gradual motion oi lo ^ L can be described analogous to a second-order-transition in a bulk system, the free energy is drawn 
to meet the branch 4crl/^ — 9>g[:^-^^L at Hsc(L) with horizontal slope. 



tension has already been stressed by Parry et al. |2lj. If (j\inc > 0, only a first order filling transition is possible (i.e., 
at Hs — Hsc the length £o jumps discontinuously from = to £o = oo, in our simplified treatment). Using Eq. Q 
in Eq. © yields 

One should not worry about the fact that for Hg Hsc this free energy excess AFs —oo, because AF in Eqs. ||2J), 
(|5|) is of order unity only, rather than scaling with any power of the linear dimension of the system. For the filling 
transition, the relevant free energy scale is AcrL^, if for Hs < Hsc we have an interface of area (2L)^ in the system. 
The free energy depression per unit area resulting from Eq.|31is of order [{Hs — Hsc)L'^]~''^ and, hence, for \Hs — Hsd of 
order the divergence in Eq. becomes problematic. Taking the limit L — s- cx) first, and then letting Hs — > Hsc 
obviously poses no problem: the free energy per unit area stays 4\/2/s (iJs , 7") for Hs < Hsc- 

We suggest now that in a large but finite bi-pyramid the behavior of the surface free energies gets modified as 
schematically shown in Fig. U] For H > Hsc{L) the interface is located at = L, in the basal plane of the bi- 
pyramid, and the free energy is reduced by a line tension contribution, / — AaL^ + &cr[^a^L. Note that, in general, 
the line tension of an interface in the basal plane of the bi-pyramid, (Jy^^^, where two planes (with surface fields +Hs 
and —Hs) meet under an angle 2a (Fig. 12 left part), can be expected to differ from the line tension aune of an 
interface that meets a flat surface under an angle a (Fig. [21 right part), with a surface field +Hs on both sides of 
the interface. It is the latter quantity, however, which determines the scale of the depression AFs below the leading 
variation AV2L^fsiHs) = iaL^ + AV2L^f,{Hs - Hsc). Since properties like the total magnetization of the bi-pyramid 
change continuously, when £q increases from small values to -^q = L, it is assumed that the transition from the state 
with broken symmetry (the interface being located at Iq or 2L ^ £q, respectively) for Hs < Hsc{L) to the symmetric 
state where £q = L la a, second order transition, implying that the two branches of the surface free energy meet at 
Hs = Hsc{L) with a common tangent. In our notation (and in Fig.[2Jl we have allowed for a shift of Hsc{L) due to 
finite size from its asymptotic value Hsc = liniL^oo -ffsc(i)- Of course, in reality we must expect that for finite L 



5 



there is a rounding of the transition in addition to the shift, and hence there does not exist any value Hsc{L) where 
singularities of the considered model (Figs. ^ and occur for finite L, but this rounding of the transition can only 
be allowed for if statistical fluctuations are taken into account. 

If one could take the result for the free energy depression through the formation of an interface, Eq. literally, 
the resulting behavior of the surface free energies would even be slightly more complicated then conjectured in Fig.|31 
In fact, if we consider the surface free energy for Hg < Hsc{L) explicitly 

f{Hs) - AaL" + aV2L\H, - H,,)fi + AF, , H, < Hsc{L) (6) 
and insert Eq. |SJ), we recognize that f{Hs) exhibits a maximum at 

Hs{L)raax ~ H sc = ^j'l,' 

which corresponds to a value — L. Using this result in Eq. ©, we readily find f{Hs) — 4:aL^ + 8criinci at the point 
marked by a dot in Fig. 13 as ff^ — > H'^^{L). On the other hand, approaching the transition from the other side ( 
Hs Hf^{L) ), we have f{Hs) — iaL'^ + 8<7[^^^L, and there is a priori no reason to assume that crii„c = CTj'jjj^, as noted 
above, because of the physical distinction between the contact line on a plane and the contact line pinned to the edge 
where the surface field changes sign (c.f. Fig. So one would predict that in the surface free energy at Hsc{L) there 
is a jump-singularity of order L. However, as will be discussed in the next section, thermal fluctuations are expected 
to smooth out this singularity and, hence, fs{Hs) is a smooth non-singular function for all L. 



B. Interfacial fluctuations at fllling transitions 

Here we recall results due to Parry et al. on the filling of rotationally symmetric (infinite) cones. These authors 
presented arguments that the dominating fluctuations of the interface are the so-called "breather modes", i.e. the 
interface moves uniformly up or down. If we denote this fluctuating height of the interface midpoint over the cone 
(or pyramid) corner as £ and interpret £o of the previous subsection as its average value, £o = {£), the probability 
distribution derived by Parry et al |2l| for the interface position £ can be rewritten as 



P{£) (X exp 



(8) 



where the correlation length describing the interfacial width due to these breather fluctuations diverges by a simple 
power law 

u « ir"" ■■''^-l- (9) 

Although Eqs. (jHJ and @ have been directly obtained for axially symmetric cones only. Parry et al. assert that 
they should hold as well for the inverted pyramid-shaped geometry considered in the right part of Fig. 13 

It is interesting to apply Eqs. ((SJ and (jSJ to such a semi-infinite inverted pyramid for the case when £q has reached 
the value £o = L. Then Eqs. © and © imply, using £o oc t~^, 

{£^) ~ {£f = (X \t\~^ (X £o (X L . (10) 

This result also means, however, that the fluctuation in the area of the interface is proportional to L as well. Since 
these fluctuations of the interfacial free energy are of the same order as the contribution of the line tension, a theory 
based on balancing surface free energy differences with the line tension alone, as sketched in the previous subsection, 
cannot be expected to be quantitatively valid. Basically, one could argue, what needs to be done is to average the 
free energy function of the previous section with the Gaussian distribution resulting from Eq. the result will then 
be a renormalized effective free energy varying smoothly with Hg, as anticipated in Fig. |31 

Eq. UlUI) yields a justification for the assumption that only the uniform "breather" mode needs to be taken into 
account while all the nonuniform interfacial fluctuations can be neglec ted. As is well known, long wavelength nonuni- 
form interfacial fluctuations can be modeled as capillary waves |33l l48l Isol Isil , and over a length scale L these 
capillary waves cause a broadening of the interfacial profile described by the following expression for the mean square 
width w'^ 

w^=wl + ^ln-, (11) 
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where wq is a (hypothetical |53) intrinsic width and B is a short wavelength cutoff of the same order as wq- The 
logarithmic variation of with L in Eq. results from integrating the mean square amplitude {\h{q)\'^) of the 
Fourier components h{q) of the deviation h{x) = i{x) — {£) of the local height of the interface £{x) from its average 
value 

mq)\')^^ (12) 

over all wave numbers q in the interval 2tt/L < q < 2Tr/B. 

The dominance of the uniform "breather mode" over the nonuniform capillary waves is not unique to the problem of 
the filling transition. Also in the problem of interface localization-delocalization transitions in thin films [sL 154. ISSl lSa. 
Is^ IssL I59I I of thickness D a related anomalous size dependence of interfacial widths was observed ^| . Specifically, it 
was found that for a fixed linear dimension L parallel to the competing walls (in the "soft mode" phase jSJ, |53 where 
the interface is unbound from the walls) the mean square fluctuation of the interface scales even quadratically with 

D 13 

w oc — > 00 , L fixed , (13) 

cr L 

corresponding to a fluctuation of the interface as a whole in the direction normal to the interface over a finite fraction 
of the film thickness. In Eq. (|13|l . is a length characterizing the exponential decay of the (short range) repulsive 
effective potential acting on the interface from the wall at z = 0, V{z) oc exp[— kz] + exp[— k(D — z)], z being the 
distance of the mean position of the interface from the left wall. In the opposite limit, a linear variation of the mean 
square fluctuation with D was found |60l l6l| 

2 2 ksTKD 

w = Wq ~\ — h const , L ^ 00, D Wo ■ (14) 

For a cubic geometry D — L we note that Eqs. H13[l. ( I14|l exhibit a smooth crossover characterized by 

oc — — kL , L ^ 00 , (15) 

(T 

which is the same type of relation as found above in the context of the filling transition, Eq. (|10l) . Unfortunately, 
when the interface average position coincides with the bi- pyramid basal plane (Fig. 12 left part), the effective interface 
potential is not known to us, and hence we cannot quantify the prefactor in the relation w'^ — {P) — {l)"^ oc L in this 
case. 



C. A phenomenological Landau-like theory for the phase transition of the Ising bi-pyramid 

As discussed in the previous sections, the phase transition sketched in Fig. cannot be understood solely from a 
macroscopic balance of surface and line free energies, but interfacial fluctuations must be taken into account, and the 
dominating fluctuation that needs to be considered, is a uniform fluctuation of the position £ of the interface around 
its average position in Fig. |21 (left part). However, to a first approximation, £q is related to the total magnetization 
per spin, m, by 

m/TOf, = [l-(£o/i)^], (16) 

where rrih is the bulk magnetization of an (infinite) Ising lattice at the same temperature. Here we use simple geomet- 
rical relations, noting that the volume of the total bi-pyramid is SL'^/S, the volume which has opposite orientation 
of the magnetization is 4£q/3, and we have neglected any surface contributions to the magnetization, assuming that 
the local magnetization is everywhere ±mf, in the bi-pyramid, right up to the surfaces and to the interface. We shall 
discuss corrections to this approximation below. 

Being interested in \m\/mh <SC 1, it makes sense to transform from £q io £q = L — io, i.e. we count the interface 
distance from the basal plane rather than the lower pyramid corner, to conclude that in this limit m/mf, w Mq/L. 
Thus we conclude that ((m) — m)/nrn, sa 3(£o ~ ^)l and hence Eq. (jSJ can be rewritten as (in the following we take 
mf, = 1 as the unit of the magnetization per spin) 



P{m) exp 



£2 



(17) 
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Remembering that ^± cx l^l"'^^^ [Eq. (O] this is equivalent to 

P{m) « exp [-consti^(-t)(r7i - (m))^] 
Comparing this expression with the general fluctuation formula [4^ l6^ 

V{m — {m)Y' 



P{m) (X exp 



(18) 



(19) 



V = SL'^/S being the volume of the system, we immediately conclude that the susceptibility per spin x at the transition 
of Fig. ^satisfies a Curie Weiss law for i < (i.e., Hs < Hgc) 



X'xL/\t\. 

Note that Eq. (|20|l holds only for L'^\t\ ^ 1, because only then Eq. (|18|l is sharply peaked at 



(20) 

i (m), and the 

Gaussian approximation for P{m) holds. It is also interesting to note that the "critical amplitude" |63| r_ in the 
power law x = ^-\t\~"' is proportional to L, i.e. divergent in the thermodynamic limit. However, this fact is trivially 
understood, since the bulk magnetic field creates a Zeeman energy H{m){8L^ scaling with volume, while shifting 
the interface between the oppositely oriented domains costs an energy proportional to the interface area (of order i^) 
only. 

For Hs > Hsc (i-C, i > 0) we expect that the susceptibility per spin, x, in the analogous relation where (m) = 0, 



P{m) (x exp 



IvmyikBTx) 



(21) 



also scales proportional to L, because the above argument with the Zeeman energy remains valid. Hence it is tempting 
to assume that there holds a Curie- Weiss law analogous to Eq. (I20() also for t > 0, and we suggest therefore that 



P(m) cx exp [-const L^|t|m^] , t > 0, L'^\t\ > 1 



(22) 



in analogy with Eq. I|18(l . Now we also remember that for t < there is a symmetry with respect to the sign of the 
magnetization, so Eq. (|18|l for H = really needs to be replaced by an expression that is symmetric with respect to 
the sign of (to) — ±too, 



P{m) cx - {exp [—const LF'{—t){m — toq)^] + exp [—const L'^{~t){m + toq)^] } 



^exp 



-const L^{-t) 



(to — TOo)^(to + TOo)^ 



4mn 



for 



Too and L^\t\ > 1. 



It then is tempting to interpret Eqs. H22|l . H23|l as limiting cases of a Landau- type theory 



Pirn) cx exp 
with an effective free energy density fL{m) [63| 
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/l(to) = /o + + 2^Lm'^ - Hm, 



where we now have added the magnetic field. Moreover Eq. I|25|l can be rewritten as 



fL(m) 



fo + — m„ 



4Ltoq 



(to — m.of'im + toq)^ 



Hm 



(23) 



(24) 



(25) 



(26) 



and thus Eq. (^^1) immediately leads to Eqs. if r = rgi, rp being a constant. Since 



-r/{LuL), we 



recover a mean-field exponent /3 = 1/2 for the power law toq = as expected for a Landau type theory. Of course, 

such a mean-field relation for toq is consistent with the jump of the magnetization expected in the thermodynamic 
limit (Fig. ^| only if the critical amplitude B diverges as L tends to infinity. 

The considerations discussed so far do not give any clue where the basic nonlinearity responsible for deviations of 
/l(to) from Gaussian behavior comes from. Disturbingly, the ^-dependence of toq (supposedly valid for toq/to;, <C 1) 
is inconsistent with that which follows when one uses the expression £o = (^ being a critical amplitude) from 
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Eq. ^ in Eq. The approach of mo/mi, to its saturation value unity goes as mo/mb — 1 oc (L|i|)~'^. Assuming 

that the term ui^rn^ describes physical effects due to the existence of corners, which take a fraction of \/V oc of 
the volume, a plausible assumption is ul = u/L'^ (= u/L''- in d dimensions). As a consequence, one predicts 




m„ = L,p(-t)V^ (27) 



Of course, the condition niQ/mi, <C 1 requires L'^\t\ ^ 1. As we shall see below, in this regime all singular behavior is 
smeared out due to finite size rounding and, hence, Eq. I|27|) is not directly observable. The same problem occurs for 
the critical isotherm, which follows from Eq. 1251) for r = 0, , ^ 0, 



m] 



\ dm 



= ULin-^ - H = u{m/Ly - H = (28) 

t,H 



as 

mo,t=oiH) = L{H/uy/^ , S = 3, (29) 

i.e. the critical amplitude of the power law mo(i?)|t=o = DH^/^ scales again proportional to L. If we generalize the 
problem to hyper-bi-pyramidic geometry in general dimensionality d, the result T± oc L remains unchanged, while the 
other critical amplitudes become B oc D oc L'^''^ [sj. Finally we note that //.(mo) = /o + rvn? /{AL) = 

fo — L{rotY / {Au) for i < 0, as expected from Fig|2| For \t\ of order the depression of fLi'm) relative to /q is 

only of order L~^, i.e. negligible on the scale of Fig.O Only for \t\ ~ 1/L a free energy of order 1/L is obtained. 

The concept that the dominant statistical fluctuations are the "breather modes", i.e. fluctuations of the uniform 
magnetization, fluctuations with a zero-dimensional phase space j2l| , is reminiscent of the behavior of Ising-like systems 
in high dimensionalities, d > 4, which exhibit mean-field critical behavior |65l | but nevertheless for a description of 
finite size rounding these variations of the uniform magnetization need to be taken into account. In brief, the statistical 
mechanics of this latter problem is formulated |43l l44l l45l | in terms of a partition function, 

Z = J dmexp[-L'^f{m)/kBT] , (30) 

assuming a d-dimensional hypercubic lattice of linear dimension L with periodic boundary conditions, where (/o is a 
constant) 

/(m) = /o + irm^ + ^um'^ _ ffm , r = r't, t = T/T^ - 1 . (31) 
With the magnetization distribution, 

PL{m) = exp [-L''f{m)/kBT] , (32) 

its moments are then calculated as 

(m'=) = j m''PL{m)dm (33) 

In this problem, however, m is a constant and does not depend on L, nor does any other L-dependence appear in 
/(m) as given in Eq. 131|l . With a little algebra [isL li^. liH it is then straightforward to show that for _ff = the 
moments (m'^) (for k even, odd moments all vanish) scale as 

{m'') = L-^'^'^Mk{tL'^''^) , (34) 

Mk being scaling functions that can be explicitly derived from Eqs. (|30|I - H33() . 

Here we follow exactly the same procedure, the only difference being that Eq. H31|) needs to be replaced by Eq. (|25|l . 
with ul = ujL? . Thus we obtain, writing \1/ = m/mQ and considering H = for simplicity 



exp 

+ 00 



PLi^) - ' "^'^ ^ (35) 
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It is seen that this distribution depends only on a single parameter, namely, 



KbTx u L/{2rQt) u 

and hence one finds that all moments are functions of this single parameter as well 

(1*1^-) ^ MtL^) , (37) 
where the scaling function is defined in terms of Eq. H35|) by simple integrals. Since 

(|mh = m^-(l*h = (^) ^ (l*l'> ^ rh.itL') , (38) 

also the scaling function rhk is a function of the scaling variable tL^ again, there is no L-dependent prefactor, unlike 
Eq. 

It is useful to consider the behavior right at < = separately, since then (if also = 0) we have simply 



PLim) = ^cxp 



~2um 



(39) 



since the L-dependence of the volume V = 8L^/3 is canceled by the L-dependence of ul, ul = u/ . From Eq. H39|) it 
is then obvious that all moments (|m|'^)t=o = /fc(0) are simple constants. If a magnetic field is included, we similarly 
conclude 

(m'') =mk{tL^,HL'') (40) 

and in particular one finds that the zero-field susceptibility at i = is finite and proportional to the volume, x(i = 
0) oc L"^. At this point we return to one - so far not really justified - key assumption of the present treatment, namely 
ul = u/L^. Equally well one could argue that the basic non-linearity, ului'^, of the effective free energy density is due 
to line tension effects rather than caused by the presence of pyramid or cone corners, and hence ul = u' /L'^ would 
result. As a consequence, itiq — ro\t\L/u', and one would still predict a diverging amplitude of the order parameter, 
mo = L^^^y/ro/u'\t\^/^, and a finite size scaling variable [Eq. ^] mlL^ /{kBTx') oc {tL'^/'^ f instead of {tL'^f . The 
free energy then would become fhi^) — fo~ (^o^)^/"', i-e. for |i|L^/^ of order unity it still would be of order L^^, as 
expected, since the regime of finite size rounding corresponds to total free energy differences in the system of order 

ksT. However, considering the distribution PL{m) for t ~ Q we would obtain PL{m) = Z^^ cup 



exp 



2"' ,Lm'^ 



gj,^^^..., , implying that for t = the moments scale as (iTil*^) cx L . As will be demonstrated in Sec. 

Ill, such a behavior clearly contradicts observation. 

It is clear that the treatment presented so far is extremely simplified, and one needs to discuss various corrections. 
For the Landau theory of phase transitions in the bulk at high dimensionalities [Eqs. H3U|) - (|34|l ] one knows that 
nonuniform order parameter fluctuations are the leading source of higher order correction terms to Eq. (|34|l [66l| . 
In our problem, the analogous fluctuations to consider would be nonuniform fluctuations of the interface. These 
fluctuations are correlated in a volume in d dimensions, and the motion of the interface back and forth in such a 
correlated volume would cause a magnetization fluctuation of the order of to^^^ . Hence we conclude that a contribution 
Ax to the susceptibihty results, 

kuT^X = - (|m|)2) cx mgel (x \t\-^l'' (41) 

We note that in d = 3 this divergence has a stronger power than the leading Curie- Weiss term, Eq. (|2U|I . but it has 
a critical amplitude which is of order unity rather than of order L. In the regime \t\l? 3> 1 the asymptotic behavior 
predicted by Eq. (|20|) should dominate, but for \t\LP' of order unity this correction may be non-negligible. Denoting 
the leading result of Eq. H20() by xo, we have in d — 3 

X = xo + Ax cx L\t\-^ + const \t\-'^l'^ = ^1*1"^! + const {\t\L^Y^I'^\ (42) 

Thus in the region of the finite size rounding of the transition, the nonuniform fluctuations yield corrections that are of 
the same order as the corrections that would result from the scaling function mfe(iL^) in Eq. (|38|l . Hence one cannot 
expect that moments (ItiI*^) calculated from Eq. H35|) are quantitatively accurate, therefore, we have not bothered to 
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work them out in full detail. This result must be expected, of course, because the result cx means that for 

\t\L'^ of order unity is of order L, the whole length scale of the Ising bi- pyramid. 
It is also is of interest to consider the generalization of Eq. H42(l to arbitrary dimensionality d, which yields 

X = xo + Ax = L\t\-'[ro + const /L\t\^''^'-'^] (43) 

This result shows that for "hyper-bi-pyramids" in d > 3 Ax indeed becomes a correction, smaller than the terms 
resulting from the scaling function in Eqs. I|35|) - (|37|l . Conversely, for d < 3, the corrections in the finite size scaling 
limit of order unity are proportional to L'^(3-d)/2 hence larger than the leading term of order tq. This 

marginal role of d = 3 may indicate possible logarithmic corrections to finite size scaling. For < 3 we expect a 
nontrivial description of finite size effects with non-mean-field exponents. Monte Carlo studies of the d = 2 square 
geometry with competing edge fields '6^ corroborate this conclusion. 

In addition to the effects of nonuniform interfacial fluctuations, there exist also corrections due to walls, edges and 
corners that invalidate the simple relation between the interface height io and the magnetization, Eq. (|16(l . E.g., 
from the regime where in Fig. El a negatively oriented domain meets positive surface fields, we expect a correction 
4V2(L2 - el){AmJmo)/{8L^/3) to Eq.UHl so that 

^^l-(V^)^-(3/V2)^(i-||) , (44) 
nib nib \L L-^ J 

AttIs being a surface magnetization difference (per unit surface area). Additional corrections (of order L~^) may 
result from the edges where the triangular surfaces of the pyramid meet. Thus a phenomenological relation between 
y = m/irib and x = (q/L is 

y = 1 - Aqx - Aix"^ ~ A2X^ , (45) 

where Aq,Ai and A2 are phenomenological constants. 
When £0 is close to L, Eq. (|44|l can be simplified as 

^^3/L_4_^Am,L_4y 
nib V ^ ™b ^ / 

which implies that in the linear relation between ((m) — m)/mb « 4(£o — ^)/^ used to justify Eq. (|17|l a correction 
term (of relative order ■\/2(Ams/mf,)/L) enters, giving rise to the replacement of the factor L^|t| in Eq. ((TH|) by a 
linear combination of terms L'^\t\ and L\t\, causing thus additional corrections to finite size scaling. 



III. SIMULATION RESULTS 



Qualitative evidence for the actual occurrence of the transition sketched in Fig. is presented in Fig. 0] showing 
two snapshots of the Ising bi-pyramid in the two "phases" caused by different values of the surface field. Note that we 
use an Ising nearest neighbor Hamiltonian 

bulk surfaces upper surfaces lower surfaces all spins 

n^-JJ2S^SJ-Js S^SJ-Hs S^ + Hs Y S,~H J2 (47) 

where the exchange constant is weakened if both spins i, j are in a surface plane, Js = J/2. One can see that for 
J/ksT = 0.45, Hg/ J = 1.00 the magnetization is still predominantly negative, as anticipated in the schematic drawing 
of Fig. 121 (left part). In the negative domain only small clusters of positively oriented spins occur, and vice versa. For 
J/ksT — 0.25 and Hs/ J — 0.77, however, there is no majority of either positive or negative spins, as far as one can 
tell this from viewing the pyramid surfaces. 

A more quantitative characterization of the transition is provided by contour diagrams (Fig. jsl . Fig. jsl shows 
that the schematic view of an interface at a height io over the pyramid corner (Figs. ^ and j^J is a very crude over- 
simplification of reality: rather the interface is fairly broad, spread out over a thickness of many lattice spacings. 
Moreover, the interface is strongly bent and not at all horizontal. One can also see that the interface is not hitting the 
external surfaces under a well-defined contact angle. Rather the midpoint contour m{x, y, z) = gradually bends over 
tangentially towards the external surfaces, see Fig.[Sli. Some crowding of contours near the external walls is always 
seen, implying that the corrections discussed in Eqs. will make a substantial contribution (62(1 . Also in the 
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FIG. 4: Snapshot pictures of the state of Ising bi-pyramids with L — iO and two surface fields, (a) Ha/ J — 1.00, J/ksT = 0.45, 
and (b) Hs/J ~ 0.77, J/ksT = 0.25. The local magnetization rrii at the pyramid surfaces is coarse-grained over each triangle 
of closest neighboring spins, thus taking on the values rrii — — 1,-1/3, 1/3, 1 with i being the considered lattice site of the 
respective surface plane. 
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FIG. 5: Contour plots presenting curves of constant magnetization (as shown in the key of each figure) as a function of position 
(z and X or y, respectively, choosing the coordinate origin in the center of mass of the bi-pyramid) for L — 26, ksT/J = 
4, Js/J = 1/2, and Hs = 0.6 (a) and 0.8 (b). 



symmetric situation (Fig.|3)), where the interface is flat and not bent, and the contour m(x,y,z) — does coincide 
with the basal plane, one can see that the effective width of the interface is quite broad. Due to this interfacial 
broadening, we expect that the details of the singular shape of the system (external surfaces with competing surface 
fields meet at z = under a sharp angle, 2a = 90*^ here) do not matter, and if the bi-pyramid edges would be rounded 
away by a smoothly curved behavior, we should still observe the same behavior as in the present study as long as the 
radius of curvature in these smoothly curved regions is a finite constant, independent of L. 

We next turn to a Monte Carlo test of the free energy constructions discussed in Fig. |21 For this purpose the 
surface free energy difference fs{Hs) is needed, and in order to find this quantity we apply thermodynamic integration 
methods 'g^, as done in our recent study of wedge filling |8(il |. 

Writing fs+{Hs), fs-{Hs) for the surface excess free energies of the bulk phases with positive (-1-) and negative (-) 
magnetization and using the symmetry relation for the Ising model fg^{Hg) = /s+(— i/s), we find that the required 
surface free energy difference can be written as fs{Hs) — fs+{Hs) — fs+{—Hs). Recalling the relation from surface 
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0.001 , 0.002 0.003 

FIG. 6: (a) Local magnetization per spin, Ms, in the surface plane [traditionally this quantity It^I is denoted as Mi, to 
avoid confusion with the surface excess contribution nis, to the total magnetization, of. Eq. (I44ll 1 . plotted versus the local 
surface magnetic field, using parameters ksT / J — 4, Js/J = 1/2, and linear dimensions L in the range 20 < L < 100, as 
indicated, (b) Surface free energy difference Aa{Hs)/ J plotted vs surface magnetic field Hs, as obtained from the data in part 
(a) via thermodynamic integration, Eq. I49II . The full straight line shows the result of an extrapolation of Fs{H)/J to the 
thermodynamic limit (see part c). Broken horizontal straight line marks the value ITll of the interfacial tension, a, of the Ising 
model at ksT/J — 4. The inset shows the extrapolation of the apparent plateau values (reached at Ha = 0.9) of Aa/J versus 

. Arrow on top shows Hsc = 0.76. (c) Finite size extrapolation of the surface free energy difference Aa{Hs)/ J plotted vs. 

, for four different choices of Hs, as indicated. Arrow shows the expected value of Aa(Hsc)/J at the phase transition. 



thermodynamics of fcrromagnets 

MMHs) = -(|^) , (48) 

where Ms^{Hs) is the local magnetization per spin in the surface plane of an Ising ferromagnet with positive magne- 
tization in the bulk, subject to surface field Hs, we recognize that the required free energy difference can be written 
as 

MHs) = 1"' Ms+{H'JdHi (49) 

Fig.Et shows a plot of Ms versus Hs for the Ising bi- pyramids with various linear dimensions L as used in the present 
study. In principle, for an accurate estimation of fs{Hs), one should use not bi-pyramid surfaces but rather surfaces 
of large parallelepipeds, (oriented with the same angle a relative to the plane 2; = as studied here), where effects due 
to edges and corners could be avoided by using suitably "staggered" periodic boundary conditions. Such a study has 
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not been attempted here, rather we work with the bi-pyramid geometry throughout, but we then have to carefully 
consider the finite size efi^ects, which are indeed quite pronounced (Fig. |S1d,c). Here Aa{Hs) — ^/2fs{Hs), to account 
for the fact that the total surface area of the triangular facets of the pyramid is \/2(2L + 1)^ in our model, while the 
total surface area of the basal plane is (2L+ 1)^, measuring lengths in units of the lattice spacing. The corresponding 
value of the interface tension of a planar interface in the Ising model at ksT/J = A, a/ J — 0.3924 JJ^, is indicated by 
a dashed horizontal line. In principle, we expect that Act/ J near H^c is a straight line, which intersects ct/J precisely 
at Hsc- In fact, the values of Aa^Hg) extrapolated to the thermodynamic limit do show such a behavior, yielding 
Hsc = 0.76. However, for all finite L the curves Act/J smoothly bend over, and reach horizontal plateaus for large 
Hg, which are substantially lower than ct/J. In terms of Fig. I^j, the occurrence of these plateaus is understood from 
the fact that for sufficiently negative Hs obeys already the symmetry Ms{Hs) = —Ms{—Hs), because the sign of 
the magnetization in the corresponding pyramid has changed when the interface has moved towards the basal plane. 
For fields Hs where this symmetry holds the integral in Eq. (|49|l yields vanishing further contributions, resulting in a 
horizontal variation of Fs{Hs) with Hg in this regime. Qualitatively, the behavior seen in Fig.|Hl3, closely resembles the 
expected behavior as hypothetically sketched in Fig.O The linear extrapolation of the saturation plateaus Act*°*/J 
versus (inset of Fig. ISJd) is nicely consistent with this picture, since the linearity of the plot asserts that the 
depression of the plateaus indeed is a line tension effect, and the extrapolated value (Act^^'/J « 0.39) agrees with 
ct/J within the statistical error. From the slope of the broken straight line in the inset in Fig. Et* we deduce the 
estimate 

8CTi'inc/^= -4.18 ±0.16 . (50) 

Unfortunately, we are not aware of any estimates of cr[i^^ for our geometry in the literature, to which our result could 
be compared. 

In the regime where Hs is small, so that the bi-pyramid has essentially a uniform magnetization, apart from the 
region close to the lower pyramid corner (cf. Fig. \^), the finite size correction to the surface free energy differ- 
ence Aa{Hs) varies proportional to L^^, as expected for a corner correction (Fig. O;). Unfortunately, a reliable 
extrapolation of Aa(Hs) in the region near Hsc would require to simulate much larger systems than was possible to 
us. 

We now turn to the description of the phase transition in terms of the moments of the distribution function 

Pl{iti) of the magnetization (Fig. [7|). Note that here and in the following m is the magnetization per spin and not 
normalized by mf,. A striking fact is the common intersection point of the {\m\) versus Hs curves (part a) for a broad 
range of choices for L, at Hsc = 0.76 ± 0.005. This value of the intersection point is fully in accord with the estimate 
resulting from the free energy intersections, obtained in Fig. ^p. The inset illustrates the fact that the slope of the 
curves {\m\) vs. Hs at Hs = Hsc increases dramatically with L; in fact, the data are roughly compatible with the 
behavior [ddml) /dHs)Hsa that one immediately derives from the scaling description, Eq. H38() . Such a rapid 

increase of the slope (d{\m\) /dHs)Hsc rather uncommon for normal second order transitions. From the curve for 
L = 40 it is already easy to guess the limiting behavior, namely (|m|) = mi, w 0.75 for Hs < Hsc, while (|to|) = 
for Hs > Hsc- Nevertheless, this jump of (|m|) resulting in the thermodynamic limit should not be mistaken for 
a standard first order transition - rather one deals here with the limiting case of a second order transition, where 
the critical amplitude diverges, and hence the critical region is exceedingly narrow. Fig. [7| demonstrates that the 
susceptibility develops a sharp peak of rapidly increasing height, as L increases. As expected from the Curie- Weiss 
law with the divergent critical amplitude, the curves do not settle at a common L-independcnt function away from 
Hg = Hsc- But the data clearly indicate a gradual growth of xl with Hs as Hsc is approached from either side of the 
transition, and the width over which this peak is rounded rapidly shrinks as L is increased. There is not a convergence 
to a delta function singularity, that would characterize a standard first order transition t iol l4l| . 

A direct analysis of these data, not requiring any bias from theory, examines the growth of the peak height with 
L, and the scaling of the peak position iJ™°^ — Hsc with L (Fig. |7|:). One nicely recognizes that also these data 
are compatible with a transition at Hsc ~ 0.76, and the relation Xmax{L) cx i"^ implies that the maximum values of 
(m^) — are of order unity, as expected on the basis of Eq. 1381) . Unlike first order transitions (where also the 

susceptibility peak height increases proportional to the volume) the width of the susceptibility peak does not shrink 
to zero for L — > oo. 

A standard method to locate critical points for various phase transitions is to check for intersections of the reduced 
fourth order cimmlant j4lll4l.l43| . 

Ul = 1- (m4)/(3(m2)2) . (51) 

Plotting hence Ul vs. Hs for various L should yield a common intersection point at Hs = Hsc- Fig. 01 shows that 
this simple recipe works here rather well again, confirming the previous estimate Hsc ~ 0.76. If one accepts Eq. H39|) 
as a description of the distribution at Hsc, one predicts for the value of U* of the cumulant at the intersection point 
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FIG. 7: (a) Plot of the absolute value (|r7i|) of the magnetization of the Ising bi-pyramid versus the surface magnetic field Hs, 
for ksT/J = 4.0, H = 0, JajJ ~ 0.5, and various linear dimensions L in the range 10 < L < 40, as indicated. Inset shows 
a log-log plot of the slope at the common intersection point vs. L. The straight line illustrates the theoretical value of the 
slope, (b) Plot of the susceptibility, calculated from magnetization fluctuations as knTxL = L^{{m^) — (|m|)^), as a function 
of the surface field Hs, including various linear dimensions L in the range 10 < L < 40, as indicated. System parameters are 
the same as in part (a). Note the logarithmic scale of the ordinate, (c) Log-log plot of the susceptibility maximum, Xmax{L), 
versus linear dimension, for the systems shown in part (b). Broken straight line illustrates the expected relation Xmax(L) oc L^. 
Inset shows the location of the maximum of xl, i?™"^, plotted vs L~^, to illustrate the convergence of //^"^ towards Hbc as 
L — > oo. (d) Cumulants Ul (Eq. I|51|l 'l. plotted vs Hs, for various L as shown in the figure, for the same system parameters as 
used in parts (a) and (b). 



the value [/* = 1 - r(5/4)r(l/4)/[3r(3/4)2] w 0.2705. The arrow in Fig. |7|i shows that this value is in very good 
agreement with the data. Interestingly, these data develop not only a common intersection point, but also a shallow 
minimum for Hs > Hsc- This is somewhat reminiscent of the behavior at thermally driven first order transitions, 
such as occur in the g-state Potts model in c? = 3 dimensions ior q > 3, where Ul is known to exhibit a very deep 
minimum f/™*" cx —L^ [t^. So the behavior of the cumulant is again indicative of a second order transition that 
is close to a first order transition. 

We now turn to a more detailed test of the finite size scaling predictions, in particular of Eq. (|38|l . Fig. |S1 shows 
scaling plots of the magnetization and the magnetization square. Using jiJ^c ^ Hs\L'^ as scaling variables in part (a), 
both branches of the scaling functions for Hs < Hsc (upper branch) and Hs > Hsc (lower branch) are combined in 
a single plot. However, the "data collapsing" on master curves in parts (a) and (c) is not really perfect, and some 
corrections to scaling are clearly seen. However, as pointed out above, {see Eqs. (|44f) - (l46|l and the accompanying 
discussion}, our scaling description has ignored contributions such as due to the surface excess magnetization iris. 
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FIG. 8: (a) Absolute value of the magnetization (|r7i|) plotted vs. the scaling variable \Hs — Hsc\L'^, using the data of Fig.|3 
and Hsc ~ 0.76 and various L as indicated in the figure. 

(b) Plot of the second moment, {m^), vs. the scaling variable {Hs — Hsc)L^ , using Hsc ~ 0.76 and including data only for 
Ha > Hac Various choices of L are shown as indicated. The straight line shows the slope of the asymptotic Curie Weiss law, 



Fig. El shows a fit of data for the magnetization to Eq. (|45|l to test the size dependence, and thus it is shown that 
indeed important corrections are present |68j . 

The simulation data in Fig. [3 were extracted from an analysis of the probability distribution Pi(m) of the magne- 
tization, m, in the finite bi- pyramid, and since Pi(m) plays a key role in our phenomenological description (Sec. II), 
we discuss Pl (to) in detail now. Fig. llOb shows that in the phase where the interface coincides with the basal plane 
of the bi-pyramid, Pl{iti) is perfectly described by the simple Gaussian, Eq. H21() . From Fig. UJi we have already seen 
that the fourth order cumulant rapidly tends to zero for H > Hsc a.s L ^ oo. The very good Gaussian fits of Fig. 
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FIG. 9: Position of the maximum of the distribution Pt(m) at Ha — 0.73 plotted vs. , for the parameters ksT/J = 4.0 
and Js/ J ~ 1/2. The bulk magnetization is — 0.750. The numerical data are fitted to Eg . 1451 and the constants Ao,Ai,A2 
are quoted in the figure. Inset shows the difference 1 — m/rrit vs. on a log-log plot. 



imply that all higher order cumulants vanish as well. Of course, this behavior is plausible due to the rapid decrease 



16 



of the fourth order term in Eq. H25|l as L ^ oo, since ul oc L~^. Thus, for {Hg — Hsc)L'^ ^ 1 the second 

moment {m?) shown in Fig.lSja already contains the full information on the distribution. 

Fig.llOb demonstrates now the smooth change of PL{m) from the single Gaussian to the double Gaussian when for 
fixed L the strength of the surface field is varied. This behavior is fully in accord with expectations for second order 
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FIG. 10: (a) Probability distribution of the magnetization PL{m) for fields Hs > Hsc, outside the critical region, using the 
parameters L = 20, ksT/J = 4, Ja/J = 1/2, and Hs = 0.8 or 1.2, respectively. Curves through the data points show fits 
to simple Gaussian, as indicated in the figure. The inset shows the function r{Hs), obtained from such fits (for L = 26) to 
demonstrate the change of sign near Hsc- The data points for Hs — 0.73 are taken from Fig. Illb below, (b) Probability 
distribution PLijn) of the magnetization, m, of an Ising bi-pyramid for L — 20, ksT/J = 4, Js/J = 0.5, and H = 0. Curves 
show various surface fields Hs/J, as indicated in the key. 



phase transitions. For a first order transition (e.g., a first order interface localization-delocalization transition has been 
recently studied both for an I sing sys tem in thin film geometry with Js/J = 1.5, Hs/J = 0.25 j?^ and for models 
of confined polymer mixtures |59l |7J| ) the corresponding distribution Pl (to) near the transition has a pronounced 
three-peak structure: two peaks have nonzero positive or negative magnetization, and the third peak occurs for m — 
(e.g. see Fig. of Ref. uM fo^' explicit example.). In contrast, here we expect near the transition a single very 
flat and broad peak (whose width should not shrink with increasing linear dimension L, as emphasized in Eq. I|39|l '). 
Fig. E] therefore examines the size dependence of PL{m) in the critical region, and part (a) shows that indeed one 
can find for each L a field H'^^{L) such that PL{m) is essentially flat near m = 0, and approximately the width of 
PL{m) stays independent of L when L increases, while H'^^^L) Hgc = 0.76 as L increases. 

Ideally, one might have expected that Eq. (|39|l should hold strictly for Hg — Hgc (be., t — 0). However, the small 
variation of H'g^{L) with L that is implied by Fig. Illb does not invalidate our phcnomcnological theory of Sec. II at 
all: as is well known 111, IM, 113) IIB IS for finite systems there is no unique "pscudocritical" point, due to the finite 
size rounding different criteria to locate a "pscudocritical" point for a finite system yield results differing from each 
other (and from the true location of the critical point) by amounts which are of the same order as the rounding, i.e. 
t (X L^^ in our case. Such an argument would imply H'^^{L) — Hgc oc L^^ here. Unfortunately, our simulation data 
are not accurate enough to check this relation (and also a larger range of linear dimensions L than what is available 
for Fig. Illb would be desirable). 

If a more rapid variation of H'gc{L) — Hgc than the second power of 1/L results, it could be attributed to corrections 
to finite size scaling, some of which were identified above [Eqs. (|44|l - (|46|) ]. Despite all shortcomings that our numerical 
results still have, we consider Fig. Illb as a highlight of the present study, since it demonstrates that in the limit L oo 
at the transition point macroscopic fluctuations occur, the magnetization varies essentially everywhere in the region 
from TO = —0.5 to to = +0.5 (in a situation where the bulk magnetization is to?, « 0.75), because the interface can 
correspondingly move freely up and down. Of course, viewing the Monte Carlo simulation as a stochastic process, such 
interface "motions" are extremely slow , an d this "critical slowing down" 75] hampers severely the statistical accuracy 
of our Monte Carlo study, as expected I3J23 ■ Note that we have applied single spin flip Monte Carlo algorithms here, 
since the Swendsen-Wang algorithm |76ll77| | or related cluster algorithms are not offering any advantage in our case, 
working for temperatures distinctly below the bulk critical temperature Tcb and in the presence of nonzero surface 
flelds. 

Fig. Illb shows then the magnetization distribution PL[m) for various L at a fixed value of Hg that is definitely 
below Hgc- One can see that with increasing values of L pronounced peaks develop with a very deep minimum in 
between; in fact, it was necessary to apply so-called "multicanonical" sampling technique (see. e.g. Refs. 0,0|) in 
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FIG. 11: (a) Probability distribution PL{m), of the magnetization, m, of Ising bi-pyramids for ksT/J = 4, Ja/J = 0.5, H — 0, 
and various choices of L with accompanying choices of Ha /J for which a flat variation of PL{m) near m — was expected 
(these choices are quoted in the figure). Full curve shows the theoretical variation from Eg. 1391 Phim) ~ exp[— am*], with 
a = 2u/(3fcsr) = (l/3)u/J 30.4. 

(b) Plot of ln[PL(m)/PL(0)]/L^ vs. m, for fcsT/J = 4, Ja/J = 0.5, = 0, at fixed Ha/ J = 0.73 and various L. The quadratic 
part near m = is described by \r]\PL(rn) / Pl{Q)\/ = O.llm"^ independent of L. 



order to be able to sample more than 15 orders of magnitude in probability with sufRcient accuracy. 

We also remark that Fig. Illb is of a very different character than the corresponding distribution for a bulk Ising 
system for T < Tct 78] : there also a deep minimum in between the peaks corresponding to the two signs of the order 
parameter occurs, but it is very flat, almost horizontal, due to configurations described as two phase coexistence (e.g., 
in a bulk Ising hypercube with periodic boundary conditions slab-like domains occur). Here the minimum of PL(m) 
near m = does not correspond to a "mixed state" of the degenerate phases (the interface being bound either to 
the top or to the bottom corner of the bi-pyramid, respectively), since no such "mixed state" can exist. Rather, the 
minimum corresponds to a uniform displacement of the interface from its stable position near one of the two corners 
to the basal plane. 

Therefore the logarithm of the distribution lnP£,(m) near to = is a simple parabola, as expected from Eqs. (|24|l . 
(|25|l . This is demonstrated clearly because we have normalized inPL{m) such that all minima coincide. Dividing out 
the predicted dependence, we see that all curves near the minimum nicely superimpose. However, it is also clear 
from Fig. Illb that a fit to the form 

ln[F^(TO)/F,(0)] 1 u 

which would be implied by Eqs. H24|l . H25|l . when = u/L^ is used, is not a good representation of the data for 
large to, and higher order terms (of order to®, to®, . . .) would be required. Of course, this is not surprising at all, since 
the saturation value of the magnetization at the considered temperature is uib = 0.75, and the distribution spans the 
range from about to « —0.6 to about +0.6, i.e. close to the saturation values of to. Of course, Eq. (|52|) is supposed 
to be valid only for ]to] <C mt,. If we ignore this problem, the "best fit" values of the data in Fig. Illb would be 
—Ar/'ikBT = 0.11, 2u/3kBT = 30. Nevertheless it is reassuring that the estimates for r resulting from the fits in 
Figs[Tn|[Ill(8r/(3fci3r) = 0.56,0,0965,-0.22, for Hs = 1.2,0.8,0.73, respectively) yield a smooth curve r{Hs), which 
has a zero close to Hgc — 0.76, though this curve is clearly not a simple straight line over the wide range of values for 
Hs that is considered here. 

Now we turn, very briefly, to the behavior in non zero bulk magnetic field H (Fig. [T^ . There are no surprises: 
we find that (to) is a function of a scaled field HL for Hg > Hgc, and of a scaled field HL^ for Hg = Hsc, as 
expected from the theory of Sec. II (cf. Eqs. (|20|l . (|40|) ). If we include in the Eq. (|39|l the magnetic field, ^^(to) cx 
exp[— aTO** + SL^mH/SksT], we can readily calculate m{H) at the critical point {Hs = Hsc,t = 0) as m(H) = 
a-V46[r(3/4)/r(l/4)]{l - (6V2)[1 - r(5/4)/(3r(3/4))] + ...}, where b = Sa^'^'^^L^H/i^ikBT). The resulting slope of 
m{H) vs. L^H/J, m{H) 0.043i^iJ/J is again in good agreement with the numerical data in the inset of Fig. 1121 
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FIG. 12: Plot of the average magnetization (m) as a function of the scaled field HL, for the parameters ksT/J — 4, Ja/J ■ 
0.5, Hs — 1.2, and three choices of L as indicated. Inset shows (m) vs. HL^ at Hsc ~ 0.76. 



1.0 



0.6 



0.4 - 



0.2 



0.0 



\ -B- 




■■\ -e- 
\ A 




\-B- 

■■■A 




\ 




\ 

o L = 20 \ 


Hsc 


□ L = 26 \ 




A L = 40 % 




Slope -1 ^ 

% 









0.0 



0.2 



0.4 



1.0 



H 



FIG. 13: Plot of the inverse distance, , of the interface from the bottom corner vs. the surface magnetic field Hg. Here 
£o is obtained from an analysis of the spatial magnetization distribution in the bi- pyramid (cf. Fig. |KJ. Three lattice sizes are 
shown, and the straight line marks the prediction £q^ oc Hsc — Hs according to Parry et al. [gJl- 



and thus provides a test that the estimate of the constant a (a ~ 30 — 34, see Fig. \llh ) is reasonable. 

Finally, Fig. E| considers the variation of the interface distance when the interface is close to a corner, and hence 
the theory of Parry et al. should straightforwardly apply (i.e., Eq. should hold). Thus we plot simply vs. 
Hg to test the resulting linear variation. Indeed one recognizes that the data are nicely compatible with the predicted 
linear variation over a reasonable range of i^^, and the extrapolated intersection point with the abscissa agrees well 
with Hsc = 0.76. Of course, we cannot expect that this relation works for i^^ > 1, since the notion of an interface 
becomes absolutely meaningless when its distance from the pyramid corner becomes of the order of a single lattice 
unit, or even less. Conversely, finite size effects set in when £o « L/2. In view of the fact, that the width of the 
critical distribution of the interface fluctuating around the basal plane of the bi-pyramid is of the order of Am w 0.5, 
see Fig. Illh. the strong finite size effect seen in Fig. El is not at all unexpected. 



19 



IV. DISCUSSION 



We start by summarizing the main findings of this study: 

i. The Ising ferromagnet in a geometry with free surfaces, where surface fields are applied such that one half of 
the surface experiences a positive surface field and the other half experiences a negative surface field (of the 
same absolute strength) has no total magnetization below the bulk critical temperature Tcb down to the filling 
transition temperature Tf{Hs), see Fig.^ This happens because the boundary conditions stabilize the interface, 
separating two domains of opposite orientation of the magnetization but taking equal volume. Below Tf{Hs), 
the interface (in the limit L — > oo) has essentially disappeared (it may be located at a finite distance of order 
unity close to either the top or the bottom corner if one specifically considers the bi-pyramid geometry as done 
in Fig. ^ . Alternatively, this transition may be driven by variation of the strengths of the surface field Hg 
through its critical value Hsc{T), note that Hg = Hsc{T) simply is the inverse function of T = Tf{Hs) in the 
(T, Hs) plane. This transition occurs in the way described here only if the parts of the surface where the surface 
field has the same sign are all adjacent to each other. E.g., for a bi-pyramid with boundary condition where the 
sign of the surface field alternates from one triangular surface to the adjacent one on the same pyramid, this 
condition would be violated, and no such phase transition occurs: rather one observes only a rounded phase 



transition near the bulk transition point 79] . 

ii. As described in Fig. 13 in the thermodynamic limit the transition can be described in analogy with bulk first 
order transitions, where an intersection of the two branches of the free energy which describe the phases occurs. 
However, here, both branches are surface free energies, scaling like the surface area (oc L^) rather than the 
volume (oc L^). Although in the limit L —> oo the transition is characterized by a discontinuous jump in the 
magnetization (cf. Figs.^andHi), it nevertheless is a second order transition, if the line tension of the boundaries 
of the interface where it meets the walls (Fig. |21 right part) is negative. This negative line tension makes it 
energetically favorable to stabilize already a domain of the minority phase for < H^c with a mesoscopic 
linear dimension £q, whereby £o in the limit L — *■ cxd diverges continuously, io oc {Hsc — Hg)^^. Thus in a 
small interval of Hg close to Hgc in Fig. [7^, or in a small interval of temperature near Tf{Hg) in Fig. ^ (right 
part), the magnetization (|to|) of the system varies continuously from its saturation value rrib to zero. In the 
thermodynamic limit, however, the width of this interval shrinks to zero. On the other hand, this width over 
which the smooth variation occurs is of the same order as the width of the interval over which the finite size 
rounding of the transition occurs, namely of order . Therefore the power law, Eq. (|27|l . predicted by a simple 
Landau-like theory that ignores finite size rounding, is nowhere clearly observable. In the finite size scaling plot 
(Fig. ISK) one cannot identify a branch with a slope 1/2 on the log- log plot on which the curves collapse. This 
is prevented by the saturation of the order parameter and, thus, the curves bend over to fiat plateaus. This 
saturation is not described by the theory we have developed here, and it clearly violates the finite size scaling, 
as is evident from Fig. 

iii. A particular interesting behavior exhibits the total susceptibility of the system. Figs.Qj andlSja imply, in accord 
with our theory, that the susceptibility x shows a Curie- Weiss- like divergence for Hg > Hgc, Eq. H20() . The 
region of bulk fields, where this divergence is observable, shrinks to zero like 1/L, because the critical amplitude 
in Eq. H2U|) varies like L. There is a remarkable asymmetry between the behavior of the susceptibility in the 
regime Hg > Hgc (where no total magnetization occurs) and the regime Hg < Hgc, however: in the latter 
regime, the total magnetization (|m|) essentially reaches its saturation value in the interval Hgc — Hg oc l/L^, 
and this is the same regime over which the Curie- Weiss-like divergence of the susceptibility is rounded off. For 
{Hgc — Hg)L^ ^ 1, however, {\m\) is almost identical to its saturation value, and the susceptibility converges 
towards the (small) susceptibility of a bulk Ising system, independent of size. The maximum value of the 
susceptibility scales like the system volume, L^, as in a first order transition (Fig. [7|:), but unlike the latter the 
shape of the susceptibility maximum does not converge to a delta- fmiction singularity (Fig.[7|D). 

iv. A very special behavior is detected for the probability distribution of the order parameter fFigs. [Tni and[TT |l . On 
a scale of \Hgc — Hg\ oc the shape of this distribution changes from a single Gaussian peak (for Hg > Hgc) 
to a double peak distribution (for Hg < Hgc), and the inverse width r{Hg) vanishes linearly with Hg as Hgc is 
approached from above (see inset of Fig. llUb). This transition from single to double peak shape happens via 
distribution Pt(m) oc exp[— am^], with a coefficient a that is independent of L (Fig. ITTk ). This broadness of 
the distribution implies that in the limit L ^ oo at Hg = Hgc fiuctuations of the magnetization occur which 
have a macroscopic, size-independent amplitude. The standard statement of statistical thermodynamics, that in 
the thermodynamic limit the relative magnitude of fluctuations (| Am|/mb) is negligibly small is not at all true 
here. Again the behavior is completely different from a standard first-order transition: at the latter, the system 
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would jump between m = and m = —rrib and +mi,, whereas at the pathological second order transition found 
here the magnetization can fluctuate over a finite fi'action of the interval between —mi, and +TOf,, characterized 
by the distribution of Fig. Illb . At a first order transition, we would instead have three delta functions (at 
m — ±mi, and m — 0, respectively) as an order parameter distribution. However, this particular behavior is 
easily accounted for by the Landau-type theory of Sec. II. In particular, for Hg < Hsc the simulations confirm 
the behavior XhPl oc L^^m? near the minimum at m = fFig. II Ih). 

V. In the regime for Hg < Hsc, when the system is large enough so that |m| is close to its saturation value m^, the 
variation of the interface distance £q counted from one of the corners is found to go as £q oc {Hsc — Hs)^^, as 
predicted by Parry et al. 0|. When £o becomes comparable to L/2, this divergence is rounded off, as expected 
from the behavior of PL{m), since then a crossover from the theory of Ref. 21] to the behavior described by our 
Landau-like theory occurs. However, the detailed behavior in this crossover regions is not yet fully understood. 

vi. We now discuss the extent to which similar behavior can be expected to be found for real systems, such as 
the liquid-gas transition in a suitable cavity, where half of the surface area has an energetic preference for the 
liquid (such that "incomplete wetting" of the liquid at the wall occurs) and the other half prefers the gas (i.e., an 
"incomplete drying" boundary condition). In this case one neither has a precise symmetry between liquid and 
gas in the bulk, nor can one expect a precise antisymmetry between the interactions at the two types of walls. 
Therefore the phase transitions will be shifted somewhat away from the chemical potential value at which phase 
coexistence can occur in the bulk. A similar smooth interpolation between "capillary condensation"-like behavior 
and "interface localization-delocalization" transitions has also been found for a model of a polymer blend in a 
thin film geometry confined between parallel plates at which surface fields act which do not have a particular 
symmetry (80i] . It is encouraging that wetting phenomena in systems where chemically distinct substrates meet 
find increasing theoretical [2^ and experimental attentionls^- An important constraint though is that for the 
fiuids one has to consider the grand-canonical ensemble where the fluid in the cavity can exchange particles with 
a reservoir, and similarly for the binary fluid in the cavity also exchanges A ^ B or vice versa must be possible, 
due to a connection with a suitable reservoir. If one considers a fluid in a cavity with a fixed total number of fluid 
particles, or a binary mixture in a cavity with flxed relative concentration, this transition of Fig.^is completely 
suppressed: this situation would correspond to an Ising system at constant total magnetization, and hence by 
construction fluctuations of the uniform magnetization then are impossible. In Fig. ^ then the conflguration 
with an interface present in the basal plane of the bi- pyramid is enforced at all temperatures. An interesting 
aspect is also the crossover between the transition studied here and the standard critical behavior near the bulk 
critical temperature (a crossover from wetting to critical adsorption 82]). 
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